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Abstract 

Complex fission phenomena are studied in a unified way. Very general reffec- 
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Oh- 



tion asymmetrical equilibrium (saddle point) nuclear shapes are obtained by solv- 



. ing an integro-differential equation without being necessary to specify a certain 



parametrization. The mass asymmetry in binary cold fission of Th and U isotopes 



D ■ is explained as the result of adding a phenomenological shell correction to the liquid 

drop model deformation energy. Applications to binary, ternary, and quaternary 
fission are outlined. 
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1 Introduction 

• rH 

X. 

There is a continuous progress in deeper understanding the large variety of fission pro- 
cesses which can be theoretically treated in a unified way pQ. Even a-decay and cluster 
radioactivities Ej can be considered among the members of this family. Particle accom- 
panied fission was discovered in 1946 (see the review j3j), but only recently, by using new 
methods of fission fragment identification based on the characteristic rotational spectra 
measured with large arrays of Germanium Compton-suppressed detectors (as GAMMA- 
SPHERE) E], the "cold" processes a and 10 Be accompanied cold fission as well as 
the double and triple fine structure in a binary and ternary fission have been discovered 
jlj . The unified approach of cold binary fission, cluster radioactivity, and a-decay was 
extended to cold ternary [H] and to multicluster fission including quaternary (two-particle 
accompanied) fission 0. In that paper we stressed the expected enhanced yield of two 
alpha accompanied fission compared to other combinations of two light particles; it was 
indeed experimentally confirmed [X0*| 111). In a cold binary fission the involved nuclei are 
neither excited nor strongly deformed, hence no neutron is evaporated from the fragments 
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or from the compound nucleus; the total kinetic energy equals the released energy. In a 
more complex than binary cold fission (ternary, quaternary, etc), neutrons could still be 
emitted from the neck, because the Q-value is positive. In this case their kinetic energy 
added to those of the fragments should exhaust the total released energy. 

The importance of scission configuration for ternary fission |T2] was repeatedly stressed 
in the past. In binary fission, it is now better understood due to a longstanding effort of 
systematic analysis (THl EH 113 EE] • 

The statical approach was widely used [13 EB1 EH! t° find the saddle point shapes within 
a liquid drop model (LDM). Usually the equilibrium nuclear shapes are obtained by 
minimizing the energy functional on a certain class of trial functions representing the 
surface equation. Such an approach shows the importance of taking into account a large 
number of deformation coordinates (it seems that 5 coordinates are frequently needed) |20] 
l2~T] . The parametrization of Legendre polynomial expansion with even order deformation 
parameters «2 n up to n = 18 was employed [d to describe various saddle point shapes 
including those very similar to two tangent spheres. 

In order to study nuclear properties we have to consider both the collective and the 
single-particle motion of nucleons. This can be done by adding a shell correction to the 
LDM deformation energy (221- Otherwise a well known asymmetrical mass distribution 
of fission fragments or the ground state deformation of the majority of nuclei could not 
be explained. 

In this paper we present results obtained with a method allowing to find a general reflection 
symmetrical or asymmetrical saddle point shape as a solution of an integro-differential 
equation without a shape parametrization apriori introduced. This equation is derived 
as a Euler-Lagrange relationship associated to the variational problem of minimizing the 
potential energy with constraints (constant volume and given deformation parameter). 
The axially-symmetrical surface shape minimizing the liquid drop energy, Eldm = E s + 
Ec, is straightforwardly obtained. Minima of the saddle point deformation energy appear 
at finite values of the mass-asymmetry parameter as soon as the shell corrections, 5E, are 
taken into account (23 121] • A phenomenological shell correction is used. Also we shall 
discuss the multicluster fission phenomena. 

2 Equilibrium shapes 

In contrast to reality, within a liquid drop model all nuclear shapes in the ground-state 
are spherical and the fission fragment mass distributions are symmetrical. Permanent 
nuclear deformations and fission fragment mass asymmetry can be explained by com- 
bining the collective (liquid drop-like) and single particle properties in the framework of 
a macroscopic-microscopic method. By using the two center shell model to describe the 
single-particle states in binary fission or the three center shell model in the ternary fission, 
one can follow the shell structure all the way from the original nucleus, over the potential 
barriers, up to the final stage of separated fragments. Particularly important points on 
a potential energy surface are those corresponding to the ground-state, saddle point and 
scission point. 
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In order to illustrate the definition of the saddle point we plotted in figure^an example of a 
LDM potential energy surface (PES) versus two deformation parameters: the elongation 
R [or the dimensionless quantity (R — Rj)/(R t — Rj) where Ri and R t are the initial 
and touching point values of the separation distance R] and the mass asymmetry 77 = 
(Ai — A2)/(Ai + A2). A statical path may be seen on this PES as a heavy line following 
the valley of the potential minimum which corresponds to 77 = 0. If we start with a large 
value of R and then follow the decreasing elongations, the bottom of that valley leads to 
increasingly higher energies up to a maximum (the saddle point) marked with a cross on 
figure ^ then the energies along the valley are decreasing until the ground state minimum 
is reached. In this example the energy E = E(R,rj) is function of two shape coordinates 
and the fission valley represents a conditional minimum of the energy [dE/drj = for 
different elongations R = Rk (k = 1,2, ...,n) and 77 = 0]. The maximum value of this 
minimum is the saddle point defined by 

OE/Ot] = OE/dR = (1) 

d 2 E a 2 E 

dR 2 dRdn ^ n (0\ 

a 2 E a^E < u \ z ) 

drjdR drj 2 

If we take another value of mass asymmetry 77 = 77^ ^ 0, for every Rj within LDM, we 
obtain a new (conditional) saddle point at higher energy (see figure El where di — d^ is 
proportional with 77), proving that the mass asymmetric distribution of fission fragments 
is not explained by a pure LDM. 



3 Integro-differential equation 

We are looking for a function p = p(z) expressing in cylindrical coordinates the nuclear 
surface equation with axial symmetry around z axis and the tips z\ and z-i- The depen- 
dence on the neutron, N, and proton, Z, numbers is contained in the surface energy of 
a spherical nucleus, E®, the fissility parameter X = Eg/(2E®), as well as in the shell 
correction of the spherical nucleus 5E°. Eq is the Coulomb energy of the spherical shape 
for which the radius is Rq = tqA 1 ^ and the mass number A = N + Z. The radius 
constant is tq = 1.2249 fm, and e 2 = 1.44 MeV-fm is the square of electron charge. The 
lengths are given in units of the radius, Rq, and the Coulomb potential at the nuclear 
surface, V s = (R /Ze)(p s , in units of Ze/R^. The surface tension and the charge density 
are denoted by o and p e respectively. The nuclear surface equation we are looking for 
should minimize the functional of potential energy of deformation 

E. + E C = 2^*5 £ P U)^?d 2 + £' - £ S£) V,dz (3) 

with two constraints: volume conservation, 

V = nRlf%\z)dz=^ (4) 
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and a given deformation parameter, 



ttR 3 

a =!lp. I ' F{z,p)p 2 dz (5) 



V oz x 

assumed to be an adiabatic variable. 

According to the calculus of variations the function p(z) minimizing the energy with 
two constraints should satisfy the Euler-Lagrange equation (see the Appendix) 

PP " - p' 2 - (Ax + \ 2 \z\ + QXV S ) p(l + p' 2 f/ 2 -1 = (6) 

or 

2aK + 3p e (j) s /b + X[ + X' 2 \z\ = (7) 
where X[ and A' 2 are Lagrange multipliers and K is the mean curvature [26J: 

K={TZ^+TZ^)/2 (8) 

with 1Z\ and 1Z 2 the principal radii of curvature given by 

U x = R rp ; TZ 2 = -R T 3 /p" ; r 2 = 1 + p' 2 (9) 

where p' = dp/dz and p" = d 2 p/dz 2 . 

It is interesting to mention that in the absence of an electric charge, the condition of 
stable equilibrium at the surface of a fluid |23 EH] is given by Laplace formula equating 
the difference of pressures with the product 2aK. 

By choosing the deformation coordinate as the distance between the centers of mass 
of the left and right fragments, a = \z c L \ + \z c R \, one can reach all intermediate stages of 
deformation from one parent nucleus to two fragments by a continuos variation of its value. 
Also a possible dynamical study, for which the center of mass treatment is very important 
|29j . may conveniently use this definition of the deformation parameter. The position of 
separation plane between fragments, z = 0, is given by the condition (dp/dz) z=0 = 0, 
which defines the median plane for a usual spherical, ellipsoidal, or "diamond" shape in 
the ground state, or the middle of the neck for an elongated reflection symmetrical shape 
on the fission path. For this choice of the function F(z,p) one has / = \z\. 

At the left hand side and right hand side tips on the symmetry axis one can write 

p( Zl ) = p{z 2 ) = (10) 

and the transversality conditions 

dp(zi) dp(z 2 ) , . 

— — = 00 ; — — = -00 (11) 

dz dz 

The equation is solved numerically by an iterative procedure checking the minimization 
of the deformation energy with a given accuracy. The phenomenological shell corrections 
to the LDM deformation energy presented in section El are used to obtain reflection asym- 
metric saddle point shapes. In fact the equation to be solved numerically is obtained from 
(JHJ) after changing the variable and function as shown below. 
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One can develop the computer code for just one of the "fragments" (for example for 
the right hand one extended from z = to z = z 2 ) and then write the result for the 
other fragment (left hand one from z = —Z\ to z = 0). For symmetrical shapes we have 
z 2 = Zp = —Z\. It is convenient to make a change of the function and variable defined by: 



u(v) = A 2 p 2 [z{v)} ; z(v) = z p - v/A 



(12) 



therefore dz/dv 



it 
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-1/A, u' = du/dv = 2A 2 p(dp/dz) dz/dv = -2App', p = y/u/A, 
Aup' 2 , l + p >2 = u' 2 /{Au) + 1, u" = d 2 u/dv 2 = d{u')/dv = -2A[dp{z{v))/dv}dp/dz - 
2Apd(p(z(v))/dz = -2A(dp jdz) {dz/dv) {dp jdz) - 2Ap{d 2 p/dz 2 ) dz/dv = 2p' 2 + 2pp" . By 
substituting into equation (JHJ) one has 



u 



1 

u 



u' 2 + 



f3XV s Ai + X 2 z p X 2 v" 
( 2A + 4A 4A2 



{Au + u l2 f 2 



(13) 



A linear function of v is introduced by adding and subtracting a + bv to 3XV S /2A. The 
quantity V s d is defined as the deviation of Coulomb potential at the nuclear surface from 
a linear function of v 

3X 

V sd = — K - a-vb 



where the constant 



is chosen to give V s d{v = 0) = 0, and 



3X 
~2A 



2A 



^V.(, = 0) 



V s (v = v p ) - a 



/v P 



where v p = Az p . Consequently one has 



u 



2- -<u 

u 
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-^A— + a +V [ b 



4A 2 



sd 



(14) 
(15) 

(16) 

{Au + u' 2 f /2 \ =0 (17) 



By equating with 1 the coefficient of v, one can establish the following link between A 
and the Lagrange multiplier A2 

A 2 = A 2 /4(6-l) (18) 
In this way u{v) is to be determined by the equation 

1 



u" - 2 - ~[u' 2 + {v-d + V sd ){Au + u' 2 f 2 } = 



a 



(19) 



where the role of a Lagrange multiplier is played by the quantity d which is taken to 
be constant instead of a. The value of the deformation coordinate a is calculated after 
obtaining a convergent solution. To the tip z = z p , at which p(z p ) = 0, corresponds v = 0, 
hence u(0) = A 2 p 2 {z p ) = 0. By multiplying with u the equation ()19j) . introducing v = 0, 
and using the relationship V s d(v = 0) = 0, it follows that u'{0) — 1/d. Consequently the 
boundary conditions for u{v) are: 



u(0) = 0, u\0) = 1/d 



(20) 
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To z = 0, at which p'(0) = (the middle of the neck for elongated shapes), corresponds 
v p = Az p and u'(v p ) = — 2Ap(0)p'(0) = 0. The point v = v p in which 

u\v pn ) = (21) 

is determined by interpolation from two consecutive values of v p leading to opposite signs 
of u'(v). The number n of changes of signs is equal to the number of necks plus one given 
in advance, e.g. for a single neck (binary fission) n = 2 and for two necks (ternary fission) 
n — 3, etc. 

Although the quantity A is not present in eq (jlfjjl we have to know it in order to obtain 
the shape function u(v). By changing the function and the variable in the eq (J3) one has 

r 3 rv pn 
A= \2Jo U ^ )dV 

and the deformation coordinate, a = z c L + z R , may also be determined by adding to 

^ f Zp 9 / x , , T r 3 f° v v - v u -dv 3 f v p. . , . , . , 
4 = 2^/ o z ,? {z)iz l V = -l^-— = —l -„)„(„)* (23) 

a similar relationship for z£. From the dependence a(d), one can obtain the inverse 
function d = d(a). 

In order to find the shape function u(v) we solve eq ()19|) with boundary conditions written 
above. One starts with given values of the constants d and n. For reflection symmetric 
shapes d^ = d R and = n^. In the first iteration one obtains the solution for a Coulomb 
potential at the nuclear surface assumed to be a linear function of v, i.e. for V s = 0. Then 
one calculates the parameters A, a, and b, which depend on the Coulomb potential and 
its deviation V s d from a linear function, and the deformation energy corresponding to 
the nuclear shape [3UIIS]. The quantity V s d determined in such a way is introduced in 
eq (|19[) and the whole procedure is repeated until the deformation energy is obtained 
with the desired accuracy. In every iteration the equation is solved numerically with the 
Runge-Kutta method. 

One can calculate for different values of deformation a (in fact for a given d^ and d^) 
the deformation energy Ed e f{ot). The particular value a s for which dEd e f{ct s )/da = 
corresponds to the extremum, i.e. the shape function describes the saddle point, and 
the unconditional extremum of the energy is the fission barrier. The other surfaces (for 
a 7^ a s ) are extrema only with condition a = constant. In this way one can compute 
the deformation energy versus di = d R . In Fig. El one can see an example of variation of 
deformation parameter and the deformation energy with di for 238 U at symmetry r] = 0. 
The saddle point corresponds to the maximum of deformation energy. 

For reflection asymmetrical shapes we need to introduce another constraint: the asymme- 
try parameter, 77, defined by 

M L -M R = Al -A 2 
1 M L + M R A x + A 2 y ' 

It should remain constant during variation of the shape function u(v). Consequently 
eq |TT)|) should be written differently for left hand side and right hand side. Now di is 
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different from d^, and so are the parameters Al and Ar. They have to fulfil matching 
conditions 

p L (z = 0) = p R (z = 0) (25) 

hence 

u L /2 (v p )/A L = uf{v p )/A R (26) 
The similar condition for derivatives p'(z) in z = 0, 

p' L (z = 0)=p' R (z = 0) = (27) 

is automatically satisfied due to eq (J2TJ) . The second derivative //'(z) can have a discon- 
tinuity in 2 = if di, ^ d R . The parameters Al and Ar are easily expressed in terms of 
rj, if we write the definition of mass asymmetry as 

2tt f l 'p 
M L = —{l + r l )=nAl i u L (v)dv (28) 
3 Jo 

M R = ^(1 - 77) = ttA^ 3 P u R (v)dv (29) 
3 Jo 

We assume that M L + M R is equal to the mass of a sphere with R = 1. It follows 

A l = (1 + t ? )" 1 / 3 A L0 (30) 

A i? = (l-r / )- 1 / 3 A R0 (31) 
where we introduced notations similar to eq f!22[) : 

A L o(i?o) = i 77 / «x,(B)(v)duf (32) 



2 7o 

The shape of a nucleus with given mass asymmetry, 77, is completely determined by the 
above written equations in which the quantities ul{v p ) and ur(v p ) are solutions of the 
eq (fTH|) . There is an almost linear dependence of r\ from the difference g?l — cZr. 



4 Mass symmetry in binary fission within LDM 

One can test the method by comparing some nuclear shapes within LDM to the standard 
results for medium and heavy nuclei. In Figure we present reflection symmetric nuclear 
shapes for binary fission of a nucleus with the fissility parameter X = 0.6 (e.g. 170 Yb), 
obtained for Ul = tir = 2 (one neck), = d R = 1.4; 1.5; 1.7, and 1.91 (for which 
a = 1.314; 1.644; 2.100 and 2.304) and a vanishing mass asymmetry 77 = 0. The saddle 
point (maximum value of the conditioned deformation energy minimum) is obtained for 
dL = 1.91, at which the shape is deformed and necked-in. 

A comparison between three nuclear shapes at the saddle point for nuclei with fissilities 
X = 0.60, 0.70, and 0.82 (corresponding to 170 Yb, 204 Pb, and 252 Cf nuclei lying on the 
line of beta-stability) is presented in Figure El One can see how the necking-in and 
the elongation are decreasing (a = 2.304; 1.822 and 1.165) when fissility increases from 
X = 0.60 to X = 0.82, in agreement with J7j. In the limit X = 1 the saddle point shape 
is spherical. The method proved its capability by reproducing the well known LDM saddle 
point shapes. 
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5 Mass asymmetry in binary fission 



Within LDM a nonzero mass asymmetry parameter (see the shapes from figure |BJ) leads to 
a deformation energy which increases with 77 as is illustrated in figureEl where 77 is replaced 
by an almost linear dependent quantity (d^ — cIr). The reflection asymmetric shapes 
plotted in figure El resulted by choosing the input parameters as follows: til = tir = 2; 
d,L = 1.40; 1.45; 1.50, and 1.60 while dft = 1.40 was kept constant, and so was X = 0.60. 
The increasing deformation energy with mass-asymmetry in figure El refers to different 
values of fissility, namely X = 0.758 for 228 Th. 

When the shell effects are taken into account a saddle point solution of the integro- 
differential equation with reflection asymmetry is obtained. In the following we shall use 
a phenomenological shell correction adapted after Myers and Swiatecki (321 • At a given 
deformation one calculates the volumes of fragments and the corresponding numbers of 
nucleons Z^a), iVj(a) (i = 1,2), proportional to the volume of each fragment. Then one 
can add for each fragment the contribution of protons and neutrons 



6E(a) = £ 5Ei(a) = Y,i 5E A a ) + SE ni (a)} 

i i 

given by 
where 

and similar eq for s(N). 



SE pi = CsiZi); SE ni = Cs(Ni 
s(Z) = Z~ 2 ^F(Z) - cZ 1 ^ 



F(n) 



/V 5 / 3 _ /v 5 / 3 

I l—l 



(n - Ni-J - n 5 / 3 + N?H 



(33) 

(34) 
(35) 

(36) 



where n G (iVi_i, Ni) is the current number of protons (Z) or neutrons (N) and A^_i, N 
are the nearest magic numbers. The parameters c = 0.2, C = 6.2 MeV were determined 
by fit to experimental masses and deformations. 

The dependence on deformation [HS] a is given by 



Lj(a) 
Ri 



(37) 



where Li(a) are the lenghths of fragments along the symmetry axis. During the de- 
formation process, the variation of separation distance between centers, a, induces the 
variation of the geometrical quantities and of the corresponding nucleon numbers. Each 
time a proton or neutron number reaches a magic value, the correction energy passes 
through a minimum, and it has a maximum at midshell. 

Results for binary cold fission of parent nuclei 226_238 Th and 23 °- 238 U are presented in 
figures El and El The minima of the saddle point energy occur at nonzero mass asymmetry 
parameters d^—dji in the range 0.04, 0.08 for the above mentioned nuclei. They correspond 
to r\ of 0.050, 0.095 which leads to Ai ~ 125 in all cases. A typical saddle point shape, for 
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232 U may be seen in figured For experimentally determined mass asymmetry p^ EIo] the 
maximum of the fission fragment mass distributions is centered on A\ = 140 in a broad 
range of mass numbers of parent nuclei. 

In order to understand correctly the figure |H1 where from the saddle point energies Egp 
of every nucleus we subtracted its minimum value we would like to give an example 

for 238 U in figure E3 In the upper part we plot the saddle point energies obtained within 
a pure LDM (see also Fig. EJ. When we add the shell corrections, the conditions of 
equilibrium are changed and in general the LDM part of the saddle point energy is not 
identical with the previous one, as may be also seen from the tables ^ and El Table 
121 shows how are changed equilibrium conditions when the shell effects are taken into 
account. There are two rows for every value of di — dp when dp — dp < 0.045: one for 
maximum value of the total E$p, and the other one for maximum value of its LDM term 

EsP-LDM- 

The minimum of the E$p is produced by the negative values of the shell corrections 
SE — 5E° which can be clearly seen in the lower part of the figure The variation of 
the saddle point energy with the mass asymmetry parameter dp — dp is almost a linear 
function of the mass asymmetry r\ for some even-mass isotopes of Th and U. The minima 
of the saddle point energy occur at nonzero mass asymmetry parameters di — dp, between 
about 0.04 and 0.085 for these nuclei. When the mass number of an isotope increases, the 
value of the mass asymmetry corresponding to the minimum of the SP energy decreases. 

As mentioned by Wilkins et al. [HB] , calculations of PES for fissioning nuclei "qualitatively 
account for an asymmetric division of mass" . From the qualitative point of view the results 
displayed in Figures |H1 and M proove the capability of the method to deal with fission mass 
and charge asymmetry. The experimentally determined mass number of the most probable 
heavy fragment [HZ! f° r U isotopes ranges from 134 to 140. The corresponding values at 
the displayed minima in Figures |H1 and El are very close to 125, which means a discrepancy 
between 6.7 % and 10.7 % for Ah- The inaccuracy in reproducing the experimental mass 
asymmetry is due to the contribution of the phenomenological shell corrections. In the 
absence of shell corrections the pure liquid drop model (LDM) reflection-symmetric saddle 
point shapes ^Zj are reproduced, and the barrier height increseas with an increased mass 
asymmetry. When the shell corrections are taken into account the LDM part behaves in 
the same manner (larger values at non-zero mass asymmetry). Only the contribution of 
shell effects can produce a minimum of the barrier height at a finite value of the mass 
asymmetry. One may hope to obtain a better agreement with experimental data by using 
a more realistic shell correction model, based for example on the recently developed two 
center shell model (SHj- 



6 Ternary Fission 

Neutron multiplicities higher than one, in the induced nuclear fission, are used to produce 
the chain reaction, on which the nuclear energetics is based. The condition of a positive 
released energy, Q > 0, in such a complex process is easily fulfilled, and the escape of 
one or several neutrons from the neck formed between the light- and heavy fragment, is 
not prevented by any Coulomb barrier. A small and narrow centrifugal barrier, due to 
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the angular momenta carried away by the neutrons, do not constitute a major obstacle. 
A charged particle has to penetrate, by quantum tunneling, a much thicker and higher 
potential barrier, leading to a long delay and to a corresponding comparable low yield. 
Nevertheless, the particle-accompanied fission (or ternary fission) was observed both in 
neutron-induced and spontaneous fission since 1946. Several such processes, in which the 
charged particle is a proton, deuteron, triton, 3_8 He, 6_11 Li, 7_14 Be, 10 ~ 17 B, 13 ~ 2 C, 15_20 N, 
15_22 0, have been detected [3]. Many other heavier isotopes of F, Ne, Na, Mg, Al, Si, P, 
S, CI, Ar, and even Ca were also mentioned. 

Different elongated shapes for ternary fissions are shown in figure ^2 For shapes with 
three fragments and two necks (n^ = = 3) when di = cIr is increased from 2.25 to 2.80 
and 7.00 the deformation a increases from 1.650 to 2.306 and 2.730. In the same time 
the elongation is initially increased from 5.234 to 5.392 and then decreased to 5.24; the 
fragment radii are 0.461/0.814/0.461, 0.592/0.753/0.592, and 0.673/0.659/0.673, leading 
to decreasing energies in units of E° s from 0.165 to 0.150 and 0.134. The configuration 
with E/E® = 0.134 is not far from a "true ternary-fission" in which the three fragments 
are almost identical: ^°Yb^|V + §V + f 4 Cr and the Q-value is 83.639 MeV. One 
may compare the above E / E° s value with the touching-point energy of these spherical 
fragments (E t — Q)/E° s = 0.239. It is larger, as expected, because of the finite neck 
of the shapes in figure ITT1 For a-accompanied fission of 170 Yb with two f^Se fragments 
Q = 87.484 MeV is larger and the touching point energy (E t — Q)/E° s = 0.103 is lower. 
A lower Q = 70.859 MeV and higher energy barrier (E t — Q)/E® = 0.147 is obtained for 
10 Be accompanied fission of 170 Yb with 33AS fission fragments. 

Figure El shows how evolves the ternary shapes for three values of the fissility (X = 0.60 
corresponding to 170 Yb, X = 0.7716 corresponding to 236 U and X = 0.8213 e.g. 252 Cf) 
when the input paramter di = d^ is increased. For smaller values of dp = d^ the 
elongation of the shapes becomes larger up to a limiting value of about 3, 4 and 4.3 
respectively; then if dp = dn is further increased it becomes slightly smaller. In what 
concerns the shapes approaching the scission into three identical fragments (figure H~3|l 
the total length increases with increasing fissility. 

Systematic calculations jSHj have shown a clear correlation between the Q-values and 
the measured yield of different isotopes for one cluster accompanied fission. For example, 
among the He isotopes with mass numbers 4, 6, and 8, 4 He leads to the maximum Q-value. 
The maximum yield was indeed experimentally observed [3] for a accompanied fission. 
Similarly, among 6 ' 8 > 10 > 12 Be, the clusters 8 Be and 10 Be give the maximum Q-values. As 8 Be 
spontaneously breaks into 2a particles it is not easy to measure 8 Be accompanied fission 
yield; consequently 10 Be has been most frequently identified. By detecting, in coincidence, 
these two alpha particles, the 8 Be accompanied fission with a larger yield compared to that 
of the 10 Be one, could be observed in the future. From 12 ,i4,i6,i8q fa e favoured is 14 C, and 
all 16 > 18 > 20 > 22 O isotopes have comparable Q-values when they are emitted in a cold binary 
fission of 252 Cf. Nevertheless, 20 O is slightly upper than the others. As a rule, if the 
Q-value is larger the barrier height is smaller, and the quantum tunneling becomes more 
probable. The stronger emission of 14 C compared to 12 C has the same explanation as for 
the 14 C radioactivity; the Q-value is larger because the heavy fragment is doubly magic. 

We should stress again that if one is interested to estimate the yield in various fission 
processes, one has to compare the potential barriers and not the Q-values. Our results 



10 



are in agreement with preceding calculations j3U] showing also preference for prolate over 
oblate shapes. Theoretically it was pointed out by Present |H] in 1941 that Uranium 
tripartition would release about 20 MeV more energy than the binary one. In spite of 
having quite large Q values [SHI, this "true ternary fission" is a rather weak process; the 
strongest phenomenon remains the a-particle-accompanied fission. 

Experiments on so-called "symmetrical tripartition" were performed e.g. using the in- 
duced fission of 235 U by thermal neutrons [32], the induced fission of 238 U by intermediate- 
energy hellium ions [12]; or spontaneous fission of 252 Cf jH], etc. An yield of 6.7 ±3.0 per 
10 6 binary fissions was reported by Rosen and Hudson [32] who employed a triple gas- 
filled ionization chamber and a suitable electronics including a triple coincidence circuit. 
Other "optimistic" results are mentioned by Iyer and Cobble j3H|, who tried radiochemi- 
cal methods of identification at intermediate energy of excitation. While at high energy 
implying bombarding heavy ions of several hundred MeV, a positive result may be 
accepted, it is not certain whether it comes out from a compund nucleus. The general 
conclusion [33] after measuring triple coincidences with detectors placed at 120° is rather 
pessimistic: except for excitation energies over 24 MeV, the true ternary yield is ex- 
tremely low: under 10~ 8 per binary fission act. By performing dynamical calculations, 
Hill arrived in his thesis and in [3E] at elongated shapes with pronounced necks looking 
more encouraging for particle-accompanied fission. It would be rewarding to perform suc- 
cessful experiments with nowadays very much improved experimental techniques, despite 
the previous rather pessimistic conclusion that "true" ternary spontaneous fission is an 
extremely rare phenomenon. 



7 Multicluster fission 

The shapes with four fragments and three necks (n^ = ur = 4) can be seen in figure ITBI 
When cLl = cLr increases from 2.30 to 2.70 and 4.00 the deformation takes the values 
2.144, 3.136, and 3.233 and the elongations are 6.077, 6.916 and 6.252. The fragment radii 
are 0.408/0.625/0.625/0.408, 0.479/0.632/0.632/0.479, and 0.608/0.616/0.616/0.608, for 
which the energies in units of E° s are 0.188, 0.216 and 0.214, respectively. The last shape, 
with E/E® = 0.214 approaches a fission into almost identical four fragments y^Yb-^^Cl 
+ fyCl + f|Ar + igAr. Again the configuration with aligned spherical fragments in touch 
is higher in energy: (E t — Q)/E° s = 0.324. Even more complex shapes can be obtained by 
further increasing the values of rii = ur. 

In 1958 it was theoretically shown (12! on the basis of the liquid drop model |H2] that for 
increasingly heavier nuclei, fission into three, then four and even five fragments becomes 
energetically more favourable than binary fission (see figure EJ)- One can take, as an ap- 
proximation of the Q value, the energy difference between the sum of Coulomb and surface 
energies for the parent (superscript °) and n identical fission fragments (superscript l ) 

n 

Q n ~(E° c + E° s )-Y:(Eb + Ei) (38) 
i=i 

where n = 2 for binary fission, n = 3 for ternary fission, n = 4 — a quaternary fission, 
n = 5 — a division into five fragments, n = 6 — a division into six fragments, etc. 
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A linear dependence of Q n on the (binary) fissility parameter, X = Eq/(2E®), of the form 

Q n /E° s ~ 1 - n 1 ^ + 2X{\ - n" 2/3 ) (39) 

has been obtained |47j . When the fissility parameter increases, fission into more than 
two equal fragments becomes energetically favored. At X > 0.426 tripartition becomes 
exothermic and for X > 0.611 the Q-value for fission into three identical fragments is 
larger than that for binary fission. The general trend, and sometimes even the absolute 
values of Q2 and Q3, are well reproduced [39J by the above equation. 

A better chance to be experimentally observed has a quaternary fission in which two 
light particles are emitted from a neck formed between two heavy fragments ED] • The 
successful experiment JTU1 \H\ on 2a-accompanied fission observed in cold neutron induced 
fission of 233 ' 235 U, confirmed our expectations. 

The possibility of a whole family of new decay modes, the multicluster accompanied fission, 
was envisaged jlHl El HIS EDI- Besides the fission into two or three fragments, a heavy or 
superheavy nucleus spontaneously breaks into four, five or six nuclei of which two are 
asymmetric or symmetric heavy fragments and the others are light clusters, e.g. a- 
particles, 10 Be, 14 C, 20 O, or combinations of them. Examples were presented for the two-, 
three- and four cluster accompanied cold fission of 252 Cf and 262 Rf, in which the emitted 
clusters are: 2a, a+ 6 He, a+ 10 Be, a+ 14 C, 3a, a+ 6 He + 10 Be, 2a+ 6 He, 2a+ 8 Be, 2a+ 14 C, 
and 4a. A comparison was made with the recently observed 252 Cf cold binary fission, 
and cold ternary (accompanied by a particle or by 10 Be cluster). The strong shell effect 
corresponding to the doubly magic heavy fragment 132 Sn is emphasized. From the analysis 
of different configurations of fragments in touch, we conclude that the most favorable 
mechanism of such a decay mode should be the cluster emission from an elongated neck 
formed between the two heavy fragments. 

In a first approximation, one can obtain an order of magnitude of the potential barrier 
height by assuming spherical shapes of all the participant nuclei. This assumption is 
realistic if the fragments are magic nuclei. For deformed fragments it leads to an over- 
estimation of the barrier. By taking into account the prolate deformations, one can get 
smaller potential barrier height, hence better condition for multicluster emission. We 
use the Yukawa-plus-exponential (Y+E) double folded model jSH E2] extended [3T] for 
different charge densities. In the decay process from one parent to several fragments, the 
nucleus deforms, reaches the touching configuration, and finally the fragments became 
completely separated. 

Within the Myers-Swiatecki's liquid drop model there is no contribution of the surface 
energy to the interaction of the separated fragments; the deformation energy has a maxi- 
mum at the touching point configuration. The proximity forces acting at small separation 
distances (within the range of strong interactions) give rise in the Y+EM to a term ex- 
pressed as folllows 



Eyij = — 4 ( — ) y/a~2i02j 



2 



4 + ~a~) 



exp(-Rij/a) 
Rij/ a 



(40) 



where 



9k = T C ° Sh (f ) " Sinh (f ) ; fk = (f ) Sinh (f 1 (41) 
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in which is the radius of the nucleus A^Z^, a = 0.68 is the diffusivity parameter, and 
Q>2i, a 2j are expressed in terms of the model constants a s , k and the nuclear composition 
parameters Jj and Ij, ci2 = a s (l — nl 2 ), cl s = 21.13 MeV, k = 2.3, / = (N — Z)/A, 
Ro = tqA 1 / 3 , ro = 1.16 fm is the radius constant, and e is the electron charge, e 2 ~ 
1.44 MeV-fm. 

In order to emphasize the strong shell effect, we have chosen in the abscisa of the figure fTTTl 
the number of neutrons of the light fragment, Nl, which in turn corresponds to a well 
defined neutron number of the heavy fragment, Njj, for a given parent and a given cluster 
accompanied fission, in our case to: Nh = 150 — Nl for 2a accompanied fission, Nh = 
148 — Nl for 3a accompanied fission, and Nh = 146 — Nl for 4a accompanied fission of 
252 Qf rjj^ ver ^j ca ] hgavy bar on each plot helps to determine the position of the magic 
number Nh = 82. In a similar way, different types of lines are drawn through the points 
belonging to the same combination of atomic numbers of the fragments Zi — Zh- The 
red full line is always reserved for the pair in which the heavy fragment is an isotope of 
Sn, (with Zh = 50 magic number). The conclusion is clear: the maximum Q- value is 
obtained when the heavy fragment is the doubly magic 132 Sn. 

In figure EH the investigated pairs are emphasized. These are j3] for the binary fission: 

102,104 Zr _150,148 Ce ^ = g^^^ 104-108 Mo _148-144 Ba ^ = Q 2 _ QQ^ 110 Ru _142 Xe = 

66), and ^Pd-^Te (Nl = 70). For cold a accompanied fission [Zj one has: jjgKr- 
^Nd (N L = 56), r l01 Sr- 5 i 2 - 147 Ce (N L = 58 - 63), i™- 104 Zr- 5 f- 144 Ba (Nl = 60 - 64), 
l° 6 - 108 Mo4 42 - 140 Xe (N L = 64-66), ^Ru^fTe (N L = 68), and ^Pd-^Sn (N L = 70). 
There is also one example of detected cold 10 Be accompanied fission of 252 Cf, namely 
lSr-i 46 Ba (N L = 58). 

On the right hand side of this figure there are plots for the new decay modes which have 
a good chance to be detected: 2a-, 3a-, and 4a accompanied fission. The corresponding 
Q-values are not smaller compared to what has been already measured, which looks very 
promising for the possibilty of detecting the 2a-, 3a-, and 4a accompanied fission decay 
modes. In fact by taking into account the mass-values of the participants, one can see 
that the Q-value for the 2a accompanied fission may be obtained by translation with 
+0.091 MeV from the Q- value of the 8 Be accompanied fission. A similar translation with 
-7.275 MeV should be made from the 12 C accompanied fission in order to obtain the Q- 
values of the 3a accompanied fission, etc. Less promising looks the combination of three 
cluster, a+ 6 He + 10 Be accompanied cold fission of 252 Cf. As mentioned above, the 2a 
accompanied fission was already observed. 

Different kinds of aligned and compact configurations of fragments in touch are shown in 
figure El On the left hand side there are three aligned fragments on the same axis, in 
the following order of the three partners: 213, 123, and 132 (or 231) and one compact 
configuration (in which every partner is in touch with all others). It is clear that the 
potential barrier for the "polar emission" (123 or 132) is much higher than that of the 
emission from the neck (213), which explains the experimentally determined low yield 
of the polar emission compared to the "equatorial" one. As it should be, the compact 
configuration posses the maximum total interaction energy, hence it has the lowest chance 
to be observed. The same is true for the quaternary fission when the two clusters are 
formed in the neck (middle part of the figure). An important conclusion can be drawn, 
by generalizing this result, namely: the multiple clusters 1,2,3, .. . should be formed in 
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a configuration of the nuclear system in which there is a relatively long neck between the 
light (n — 1) and heavy (n) fragment. Such shapes with long necks in fission have been 
considered jlH] as early as 1958. For the "true" ternary fission, in two 84 As plus 84 Ge, 
E t = 98 MeV! Despite the larger Q- value (266 MeV), the very large barrier height explains 
why this split has a low chance to be observed. 

On the right-hand side of figure El we ignore the aligned configurations in which the 
heavy fragments are not lying at the two ends of the chain. By arranging in six different 
manners the a-particle, 6 He, and 10 Be clusters between the two heavy fragments from the 
cold fission of 252 Cf, the difference in energy is relatively small. Nevertheless, the 43125 
configuration seems to give the lowest barrier height. 

The energies of the optimum configuration of fragments in touch, for the 2a-, 3a-, and 
4a accompanied cold fission of 252 Cf are not much higher than what has been already 
measured. When the parent nucleus is heavier, the multicluster emission is stronger as we 
observed by performing calculations for nuclei like 252,254 Es, 255 > 256 Fm, 258 ' 260 Md, 254 > 256 No, 
262 Lr, 261 < 262 Rf, etc. 

While the minimum energy of the most favorable aligned configuration of fragments in 
touch, when at least one cluster is not an alpha particle, becomes higher and higher with 
increasing complexity of the partners, the same quantity for multi alphas remains favor- 
able. In conclusion, we suggested since 1998 experimental searches for the multicluster 2a 
accompanied fission, for 8 Be-, 14 C- and 20 O accompanied fission. Also, the contribution of 
the single- and multi-neutron accompanied cold fission mechanism to the prompt neutron 
emission has to be determined. 



8 Conclusions 

The method of finding the most general axially-symmetric shape at the saddle point 
without introducing apriori a parametrization (inherently limited due to the finite number 
of deformation coordinates), by solving an integro-differential equation was tested for 
binary, ternary, and quaternary fission processes within a pure liquid drop model. The 
well known LDM saddle point shapes are well reproduced. The method proved its practical 
capability in what concerns fission into two, three, or four identical fragments, for which 
fission barriers given by shapes with rounded necks are, as expected, lower than those of 
aligned spherical fragments in touch. 

Nevertheless, in the absence of any shell corection it is not possible to reproduce the 
experimental data, or to give results for particle-accompanied fission. By adding (phe- 
nomenological) shell corrections we succeded to obtain minima at a finite value of mass 
asymmetry for the binary fission of 226_238 Th and 23 °- 238 U nuclei. 

Fission barriers for ternary and quaternary fission into identical fragments are lower than 
for aligned spherical fragments in touch. Our expectations concerning the possibilty to 
detect quaternary fission as 2a-accompanied fission were experimentally confirmed. 
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APPENDIX: Euler - Lagrange equation 



The functional to be minimized is the surface and Coulomb energy (given in equation 
(JHJ)) with the constraints (eqs. IHE}- We denote with Fx, F 2 , F 3 , F 4 , the corresponding 
integrands one needs to write the Euler-Lagrange equation: 



<jF x = ap(l + p' 2 ) 



1/2 



R p 



s r Rop 

^2 = 



5 

F 3 = p 2 
F - 



(p 2 ~ zpp') 



P Z F 



The derivatives are easily obtained 



d dFx 



d 



dz dp' dz 



dF\ 
dp 

PP' 



f2 



1/2 



;i+p /2)l/2 



P' 2 + PP" 



pp' 2 P " 



(l + p '2)l/2 (l + p '2)3/2 



(42) 

(43) 

(44) 
(45) 

(46) 
(47) 
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dz dp' dz 
dp 



^ = w (4 

zp) = -p- zp' (49) 



d dF 2 d , 



2p (50) 



dF A ( pdF\ , . 

V=n F+ M (51) 

Consequently the Euler-Lagrange equation can be written as 
(dF 1 d dFA R oPe <f> s (dF 2 ddF 2 \ 

° [-dp- ~ TzW ) + ~5~ \W ~ Tz W) + p ( 1 + 2/) = (52) 

leading to 

pp" - p' 2 - (Ax + A 2 |z| + 6XK) p(l + p /2 ) 3/2 -1 = (53) 

if we choose F = \z\ (hence / = \z\) and express 3Rop e /(5a) as 6X because the Coulomb 
and surface energy of a spherical nucleus within LDM are given by Eq = (3Z 2 e 2 )/(5i?o) 
and E^ = 47ri?Q<7, respectively. 

Alternatively one can obtain from this equation the equivalent relationship 

2aK + 3p e s /5 + A; + X' 2 f = ; f(z, p) = F(z, p) + f ^' p) (54) 
in which the definition of the mean curvature, K, from section El was used. 
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Figure 1: An example of a potential energy surface E^ej = Ed e f(R,r]) calculated within 
liquid drop model, with a statical path marked with the heavy line which has a cross at 
the saddle point. 
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Figure 2: The saddle point deformation energy E$p vs. mass-asymmetry within a pure 
liquid drop model. Example of 228 Th. 
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Figure 3: Deformation and potential energy vs. input parameter of the integro-differential 
equation for binary fission of 238 U within a pure LDM . 
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Figure 4: Nuclear shapes during binary fission of a nucleus with fissility X = 0.60 for 
d L = d R = 1.40,1.50,1.70,1.91. 
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X = 0.82 




Figure 5: Sadie-point shapes during binary fission of nuclei with fissility X 
0.60,0.70,0.82. 
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Figure 6: Mass asymmetric shapes during binary fission of a nucleus with fissility X = 0.60 
for 4^ dji. Pure LDM calculations. 
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Figure 7: Mass asymmetric saddle point shape of 232 U. Shell effects taken into account. 
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Figure 8: Difference between the saddle point deformation energy E$p and its minimum 
value E™™ vs mass asymmetry parameter [related to (d^ — dp)} for binary fission of Th 
isotopes. One can see the minima produced by the shell effects. 
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Figure 9: Difference between the saddle point deformation energy E$p and its minimum 
value E™™ vs mass asymmetry parameter [related to {di — dp)] for binary fission of U 
isotopes in the presence of shell corrections. 
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Figure 10: The saddle point deformation energy E$p of 238 U within a pure liquid drop 
model (top). The minimum of the E$p is produced by the negative values of the shell 
corrections (SE — 5E°)sp (bottom). 
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Figure 11: Shapes obtained by solving an integro-differential equation for til = tlr = 3, 
d L = d R = 2.25, 2.80, and 7.00. The binary fissility X = 0.60 corresponds to 170 Yb. 




Figure 12: Evolution of ternary shapes (n^ = = 3) when the input paramter di = 
is increased as shown for 170 Yb (top), 236 U (middle), and 252 Cf (bottom). The total length 
on the symmetry axis increases with increasing d^ = cLr in the upper part of the panel 
and decreases in the lower part. 
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Figure 13: Ternary shapes (til = Ur = 3) approaching a scission into three fragments 
with identical radii 170 Yb (bottom), 236 U (middle), and 252 Cf (top). 




X (fissility parameter) 

Figure 14: Approximation of Q-values for fission into equally sized fagments. 
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Figure 15: Nuclear shapes during quaternary fission of a nucleus with fissility X = 0.60 
for d L = d R = 2.30,2.70,4.00. 
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Figure 16: Q- values for the cold fission of 252 Cf vs. neutron number of the light fragment. 
The experimentally determined cold fission, a-, and 10 Be accompanied fission on the left 
hand side are emphasized. The vertical heavy bar on each graph corresponds to a magic 
neutron number of the heavy fragment Nh = 82. 
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Figure 17: Aligned and compact configurations for a accompanied cold fission (left hand 
side) and for a + 10 Be accompanied cold fission (middle). On the right hand side there 
are aligned configurations with three clusters between the light and heavy fragment for 
a + 6 He + 10 Be accompanied cold fission. The parent is 252 Cf, and the heavy fragment 
is the doubly magic 132 Sn. The corresponding energies are shown. 
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Table 1: Conditional saddle point elongation parameter dp_$p and energy E$p for a 
given asymmetry dp — dp within a pure liquid drop model for 238 U parent nucleus. 



dR-sp 


d L ~ d R 


E SP (MeV) 


1.405 


0.000 


5.822 


1.400 


0.005 


5.840 


1.400 


0.010 


5.889 


1.395 


0.015 


5.972 


1.410 


0.020 


6.029 


1.395 


0.025 


6.237 


1.390 


0.030 


6.420 


1.390 


0.035 


6.637 


1.390 


0.040 


6.888 


1.390 


0.045 


7.172 


1.390 


0.050 


7.492 


1.390 


0.060 


8.232 


1.390 


0.070 


9.115 


1.395 


0.080 


10.144 



Table 2: Conditional saddle point elongation parameter dp and energy E$p for a given 
asymmetry dp — dp within liquid drop model plus phenomenological shell correction for 
238 U parent nucleus. The shell correction for the spherical parent SE° = —5.385 MeV. 
For dp — dp < 0.045 the dp at which E$p is maximum differs from dp at which its LDM 
part is maximum. 



d R 


dp — dp 


EsP-LDM 


5E SP - 5E" 


Esp 






MeV 


MeV 


MeV 


1.385 


0.000 


5.728 


0.462 


6.189 


1.390 


0.000 


5.775 


0.413 


6.188 


1.385 


0.005 


5.769 


0.410 


6.179 


1.410 


0.005 


5.823 


0.286 


6.109 


1.395 


0.010 


5.804 


0.345 


6.150 


1.405 


0.010 


5.859 


0.250 


6.109 


1.385 


0.015 


5.936 


0.146 


6.082 


1.390 


0.015 


5.841 


0.282 


6.123 


1.385 


0.020 


6.062 


-0.064 


5.998 


1.390 


0.020 


5.953 


0.129 


6.082 


1.385 


0.025 


6.219 


-0.327 


5.892 


1.400 


0.025 


6.108 


-0.161 


5.948 


1.385 


0.030 


6.408 


-0.642 


5.766 


1.395 


0.030 


6.233 


-0.375 


5.858 


1.385 


0.035 


6.628 


-1.010 


5.618 


1.395 


0.035 


6.421 


-0.665 


5.755 


1.385 


0.040 


6.881 


-1.431 


5.450 


1.395 


0.040 


6.634 


-0.996 


5.639 


1.415 


0.045 


7.044 


-1.201 


5.844 


1.405 


0.050 


7.445 


-0.289 


7.155 


1.400 


0.055 


7.825 


0.541 


8.366 


1.395 


0.060 


8.229 


1.296 


9.525 


1.390 


0.065 


8.657 


2.403 


11.059 
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